Single-Cycle High-Intensity Electromagnetic Pulse Generation in the 
Interaction of a Plasma Wakefield with Nonlinear Coherent Structures 
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The interaction of coherent nonlinear structures (such as sub-cycle solitons, electron vortices and 
wake Langmuir waves) with a strong wake wave in a collisionless plasma can be exploited in order 
to produce ultra-short electromagnetic pulses. The electromagnetic field of a coherent nonlinear 
structure is partially reflected by the electron density modulations of the incident wake wave and a 
single-cycle high-intensity electromagnetic pulse is formed. Due to the Doppler effect the length of 
this pulse is much shorter than that of the coherent nonlinear structure. This process is illustrated 
with two-dimensional Particle-in-Cell simulations. The considered laser-plasma interaction regimes 
can be achieved in present day experiments and can be used for plasma diagnostics. 
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I. INTRODUCTION 



Present day laser technology allows us to generate ul- 
traintense laser pulses with intensities approaching 10 22 
W/cm 2 1]. The quiver energy of the electrons in the 
electromagnetic (e. m.) fields of such intense pulses is 
equal to, or greater than, their rest energy. This regime is 
achieved for l^m wavelength pulses when the pulse inten- 
sity exceeds 10 18 W/cm 2 . In this regime the relativistic 
dynamics of the electrons in the plasma inside which the 
pulse propagates, introduces a new type of nonlinear phe- 
nomena (see e. g., the review articles pj and @] and the 
literature quoted therein) that arise from the nonlinearity 
of the Lorentz force and of the relationship between parti- 
cle momentum and velocity and, at very large intensities, 
from nonlinear quantum electrodynamics effects such as 
electron-positron pair creation |3J. It was soon realized 
that such nonlinear processes can be harnessed ("rela- 
tivistic engineering" , as introduced in Ref. Q ) in order 
to concentrate the e. m. radiation in space and in time 
and produce e. m. pulses of unprecedented high intensity 
or short duration that can be used to explore ultra-high 
energy density effects in plasmas. The new possibilities 
that are made available by nonlinear relativistic optics 
in plasmas were emphasized by the results presented in 
0, where it is shown that synchronized attosecond e. m. 
pulses and attosecond electron bunches can be produced 
during the interaction of tightly focused, ultrashort laser 
pulses with overdense plasmas. The property of non- 
linear systems to respond anharmonically to a periodic 
driving force was exploited in Ref. |(J where the propa- 
gation of a high intensity short laser pulse in a thin wall 
hollow channel was shown to produce a coherent ultra- 
short pulse with very short wavelength that propagates 
outwards through the channel walls. 



Recently, a different method of generating ultra-short 
e. m. pulses was proposed in Ref. 0]. This method 
uses the interaction between a relativistic electromag- 
netic sub-cycle soliton and the density modulations of 
a Langmuir wakefield in a plasma. The mechanism en- 
visaged is based on the results of Ref. Q| , where it was 
shown, that when a laser pulse interacts with a break- 
ing wake plasma wave, part of the pulse is reflected in 
the form of a highly compressed and focused e. m. pulse 
with an up-shifted carrier frequency due to the Doppler 
effect. The pulse enhancement of the pulse intensity and 
the pulse compression arise because the electron density 
modulations in the wake wave act as parabolic relativistic 
mirrors. In the approach introduced in Ref. [7] the role 
of laser pulse is taken by a sub-cycle soliton produced by 
another laser pulse in the plasma. 

Relativistic e. m. sub-cycle solitons are formed dur- 
ing the interaction of a high intensity laser pulse with 
a plasma and a significant fraction of the pulse energy 
can be trapped in these structures, up to of 30% - 40%, 
as was shown in Refs. 0. This trapping occurs be- 
cause, as the laser pulse propagates in the plasma, it loses 
part of its energy. Since the number of photons in the 
pulse is approximately conserved, the loss of pulse energy 
leads to the down-shift of the pulse frequency 0, |t| be- 
low the electron plasma frequency (Langmuir frequency) 
ujp e = (47rn e e 2 /TO e ) 1 / 2 , where n e is the electron density, 
e and m e - electron charge and mass. As a result, part of 
the pulse energy becomes trapped inside electron density 
cavities in the form of low-frequency radiation. The typi- 
cal size of these solitons is of the order of the collisionless 
electron skin depth d e = c/uo pe . The e. m. fields inside 
the solitons consist of synchronously oscillating electric 
and magnetic fields plus a steady electrostatic field which 
arises from charge separation as electrons are pushed out- 
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ward by the ponderomotive force of the oscillating fields, 
[ToL ITU . The development of the analytical theory of in- 
tense e. m. solitons in collisionless plasma can be traced 

in Refs. [1 II El El El G3. ° na bng tim6 SCal<3 ' 
the effects of the ion motion become important and the 
ponderomotive force forms quasi-neutral cavities in the 
plasma density, named post-solitons in Ref. Ha|. Post- 
solitons were observed experimentally in Ref. |l4| with 
the use of the proton imaging technique EH- However, 
for simplicity, in the present paper we shall consider con- 
ditions when the effects of the ion motion can be ne- 
glected. 

On the basis of one-dimensional (ID) analytical calcu- 
lations, in Ref. 7] it was shown that the e. m. field of 
the sub-cycle soliton is partially reflected by the electron 
density modulations of the wake wave and that the fre- 
quency of the reflected pulse is up-shifted by a factor 27^ 
as compared to the soliton frequency and its intensity is 
proportional to 7^, where j p h is the Lorentz factor cor- 
responding to the phase velocity v p h of the wake wave. It 
was thus proven that attosecond pulses can in principle 
be generated by exploiting the soliton-wakefield interac- 
tion. 

The aim of the present paper is twofold: first, we con- 
firm and qualify the analytical results of Ref. [3 on the 
basis of two-dimensional (2D) particle- in-cell (PIC) simu- 
lations of the interaction of a breaking wake plasma wave 
with a sub-cycle soliton; second, we extend the analyti- 
cal results of Ref. to other types of coherent nonlinear 
structures - an electron vortex and a wake field, and these 
two cases are also illustrated with 2D PIC simulations. 
The parameters of the simulations presented here are de- 
liberately chosen so as to reveal unambiguously the effect 
of the reflection and the frequency up-shift, thus pro- 
viding a proof-of-principle numerical experiment. Nev- 
ertheless, one can see that these effects should occur in 
present-day experiments or can be deliberately realized 
with contemporary terawatt laser systems. 

The paper is organized as follows. The interaction of 
a wake wave with a soliton is discussed in Sec. [H] In the 
first part of this section, 111 Al we briefly recall the main 
properties of the reflected pulse, as obtained analytically 
for a soliton in the ID approximation. The results of 
the 2D numerical simulations are presented in the sec- 
ond part, 111 HI together with the comparison with the 
ID analytic results. In the two following sections IlIII and 
IIVI we present two other possible mechanisms of ultra- 
short pulse generation which are due to the interaction 
of a breaking wake plasma wave with either an electron 
vortex or with another wake plasma wave oriented in the 
perpendicular direction. Each of these sections is subdi- 
vided into two parts, the first contains the ID analytical 
model and the second - the results of 2D PIC simulations. 
Finally, in Sec. [V] the main results and conclusions are 
listed. 



II. REFLECTION OF THE 
ELECTROMAGNETIC FIELD OF A SOLITON 

A. One-dimensional theory 

In this section we recall the ID results on the form of 
the e. m. pulse reflected in the soliton-wake field inter- 
action obtained analytically in Ref. ■ 

When an intense short laser pulse interacts with a 
plasma, it induces a wakefield fl6| consisting of nonlinear 
Langmuir waves with a phase velocity v p h = (3 P h,c equal 
to the group velocity of the laser pulse. If the laser pulse 
propagates in a low density plasma, the latter is close to 
the speed of light in vacuum. The nonlinearity of the 
strong wakefield leads to the steepening of its profile and 
to the formation of sharply localized maxima (spikes) in 
the electron density E3- At wavebreak (see Ref. 
and references therein) the electron density in the spikes 
tends to infinity, 
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but remains integrable. Here X — x — v p ht, n^o is the ion 
density (equal to the unperturbed electron density n e o), 
and the ratio j3 e between the speed of the electrons and 
the speed of light varies —(3 p h to (3 p h- Close to the wave 
breaking conditions, we can write 



n e (X)~n e0 [l + \ p 6(X)}/2, 



(2) 



where A p is the plasma wavelength in the wave-breaking 
regime, S(X) is the Dirac delta-function. This density 
spike partially reflects a counterpropagating e. m. wave, 
as shown in Refs. [1 Q . 

In Ref. the reflection of a ID, circularly polarized 
soliton was considered. The soliton was described by the 
dimensionless vector potential e A/ (m e c 2 ) = A(x,t) (e y + 
ie z ) (see E3), 



A(x,t) 



2e(tts) cosh [e(ris)w pe x/c] e 
cosh 2 [e{Vls)uJ P ex/c\ — s 2 (fls) 



(3) 



where tts < w pe is the soliton frequency, and e(Os) = 

(l-0|/^ e ) 1/2 . 

The properties of the reflected pulse were derived by 
performing a Lorentz transformation to the reference 
frame where the wake plasma wave is at rest. In this 
frame the reflection coefficient 



q — ilu 



(4) 



where q = 2u) pe (2~/ p i l ) 1 / 2 , derived in Ref. was used for 
the frequency components obtained by Fourier expand- 
ing the soliton amplitude. The form and amplitude of the 
reflected pulse in the moving frame were then obtained 
by an inverse Fourier transform followed by the inverse 
Lorentz transformation of the vector potential. The ex- 
plicit form of the reflected pulse is presented in Ref. 
by Eqs. (22) to (28) and by Figs. 2 and 3. 
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The main conclusion of this analysis is that the ampli- 
tudes of the electric and magnetic fields in the reflected 
pulse are increased by the factor 7^ 2 , i. e., the pulse in- 
tensity is proportional to 7^, while its frequency is up- 
shifted by 2~f 2 h . This scaling indicates that in a tenuous 
plasma the frequency up-shift of the reflected pulse, and 
its related compression, would be so large that it could 
lead to the generation of attosecond pulses. The Lorentz 
factor 7 P h of the wakefield generated by a laser pulse in 
plasma is of the order of j p h ~ uJd/uj pe , where uid is the 
frequency of the laser pulse (driver) that generates the 
wake plasma wave (see Ref. 16] ). The frequency up-shift 
is 2^ h Vt s ~ 2^ 2 h u pe w 2j ph L0 d - Thus, for a 1 fim wave- 
length laser pulse, corresponding to the critical plasma 
density n cr = m e tj^/47re 2 w 10 21 cra -3 , the factor 2 r y p f l 
that would be required to generate an attosecond re- 
flected pulse is of order 10 3 , i. e. the density of the plasma 
must be order of 4 x 10 15 cto~ 3 . Denoting by as the di- 
mensionless amplitude of the soliton (defined in terms of 
the soliton frequency in the laboratory frame), the in- 
tensity of the e. m. field in the soliton is given by Is ~ 
{ashph) 2 x 10 18 W/cm 2 . Then, the intensity of the re- 
flected e. m. pulse is I re f w "l ph Is = a-slph x 10 18 W/an 2 . 
According to Ref. Q the paraboloidal shape of the break- 
ing wake plasma wave focuses the reflected e. m. pulse, 
further increasing its amplitude. In the case of the soli- 
ton, the enhanced scaling of the reflected wave intensity, 
Iref ~ Iph^s, which leads to 



a hph x W 18 W/cm 2 . 



B. Two-dimensional PIC simulation 

In Ref. Q the generation of single cycle electromag- 
netic pulses during the interaction of the wake Langmuir 
wave with a relativistic soliton (and in the present pa- 
per with other nonlinear coherent structures) is predicted 
with the help of a one-dimensional model. In order to 
take into account the effects of multi-dimensional geom- 
etry and strongly nonlinear plasma dynamics, as well as 
the influence of kinetic effects, we performed 2D PIC sim- 
ulations using the code REMP based on the PIC method 
and "density decomposition scheme" [l8| . 

In the simulations presented here, the grid mesh size 
is Ad/20; space and time unit is Xd and 2n /uJd, respec- 
tively. Here Xd and Wd is the driver laser wavelength and 
frequency, respectively. In the figures, the electric and 
magnetic field components are normalized to m e tOdc/e 
and the electron density is normalized to the critical den- 
sity n cr = m e Wj/ Aire 2 . 

The ions are assumed to form an immobile neutralizing 
background and thus only the electron motion is taken 
into account. This approximation is applicable because 
the typical interaction period is much shorter than the 
ion response time, e. g., in a hydrogen plasma. In the 
simulations, the boundary conditions are absorbing for 
the e. m. field and the quasi-particles. The absorbing 
condition for the e. m. fields is implemented using the 
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FIG. 1: The electric field components E z (a) and E y (b) and 
the electron density n e (c) before the interaction between the 
wake wave and the soliton, t = 61 X 2-k/uj. The soliton is 
generated by the auxiliary laser pulse which has already left 
the simulation box. The dashed circle in (b), (c) denotes the 
soliton location. 



scheme [Tjj at the cost of an additional (absorbing) edge 
in the simulation box. 

The interaction of a wake wave with a soliton is sim- 
ulated in a box with size 60 x 40A<j, including the ab- 
sorbing edges of thickness 3A<j. The results are shown in 
Figs. 11131 A single relativistic e. m. sub-cycle soliton 
is generated by an auxiliary laser pulse with wavelength 
A Q = 2Xd and dimensionless amplitude a a = 0.5, cor- 
responding to the peak intensity a 2 x I\, where 1\ — 
1.37 x 10 18 W/cm 2 x (1 fim/X 2 d ). The pulse is Gaussian 
with FWHM size (length x waist) 4Xd x 6X d - The auxil- 
iary laser pulse is linearly polarized with its electric field 
along the z-axis; it is generated at the bottom bound- 
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FIG. 2: Interaction of the wake wave with the soliton. The 
electric field component E z (a,c) and electron density n e (b,d) 
at t = 76 (a,b) and 78 x 2n/u} (c,d). Dashed line £: see next 
figure. 
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FIG. 3: Reflection of the soliton electromagnetic field by the 
electron density cusps. The electric field component E z and 
electron density n e at t — 76, 78 x 2n/u> along the dashed line 
£ on the previous figure. The first and second reflected pulses 
are marked with "I" and "2" , respectively. 



ary at t = and propagates along the y-axis at x = 20. 
The plasma wakefield, which interacts with the soliton, is 
formed by a Gaussian laser pulse, the driver pulse, with 
amplitude ad = 1-5 and FWHM size 2Xd x 12Ad, starting 
at time t = 45 from the left boundary and propagating 
along the s-axis. The driver laser pulse is linearly polar- 
ized, its electric field is directed along the y-axis. The 
plasma slab occupies the region 5<:c<35,5<y<35; 
it is homogeneous in the direction of the y-axis and it has 
convex parabolic slopes along the ir-axis from x = 5 to 11 
and from 29 to 35. This plasma-vacuum interface profile 
is chosen so as to make the laser pulse entrance into the 
plasma smoother and to avoid a fast wake wave breaking 
that could happen in the case of a sharp plasma bound- 
ary. The electron density at the center of the plasma 
slab is n e — 0.09n cr , corresponding to the Langmuir 
frequency uj pe = 0.3. The number of quasi-particles is 
3.24 x 10 6 . 

The phase velocity of the wakefield when it starts to 
interact with the soliton is v p h ~ 0.925, corresponding 
to the Lorentz factor jph « 2.63. The Lorentz factor 
7 p ft is substantially smaller than the ratio between the 
driver laser frequency and the plasma frequency because, 
in the case of short pulses, the laser pulse group velocity 
strongly depends on the pulse size |20| . 

Fig. njshows a portion of the simulation box shortly be- 
fore the interaction. The auxiliary laser pulse has already 
gone through: in its wake we see a single s-polarized rel- 
ativistic e. m. sub-cycle soliton and remnants of a broken 
wakefield at the bottom of the window. The soliton fre- 
quency is well below the unperturbed plasma frequency, 
Qs ~ 0.25uJd < ojpe- The soliton appears as a region 
of low electron density, as consistent with the fact that 
the electrons are pushed outwards by the ponderomo- 
tive force of the oscillating e. m. fields inside the soliton. 
Since the driver and the auxiliary laser pulses have differ- 
ent polarizations and the soliton inherits its polarization 
from the auxiliary laser, it is easy to distinguish the e. m. 
field reflected from the soliton in the distribution of the 
E z component. The driver laser pulse induces a strong 
wakefield which is seen in the electron density distribu- 
tion as a series of wide regions of rarefaction alternating 
with thin horseshoe-shaped regions of compression. This 
is the typical pattern of the wake of a Gaussian intense 
short laser pulse. Regions of compression correspond to 
spikes (cusps) in the longitudinal profile of the electron 
density, which is well explained by the one-dimensional 
theory |17| . These density cusps play the role of semi- 
transparent mirrors, moving with relativistic velocity. 

Figs. [21and[3]show the interaction of the density cusps 
in the wake of the driver pulse with the soliton. In Fig. [5] 
the z-component of electric field and the electron density 
are shown. The wake wave of the driver is close to the 
wave breaking regime. Each electron density maximum 
(each cusp) in the wake acts as a fast moving semitrans- 
parent parabolic mirror that partially reflects the e. m. 
fields of the soliton as it propagates through the soliton. 
The process is repeated when the subsequent cusps of the 
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electron density propagate through the soliton. Thus a 
set of short e. m. pulses is formed. We note that indi- 
vidual electrons perform an oscillatory motion, while the 
electron density cusps exhibit a progressive motion. Ac- 
cording to Maxwell equations this electric charge density, 
and the associated electric current density, determine the 
e. m. held evolution and reflection. Even though the elec- 
tron density cusp is substantially distorted as it moves 
through the soliton, it recovers after leaving the soliton. 
Surprisingly, this transient distortion of the cusp when 
crossing the soliton does not prevent the formation of 
well pronounced single-cycle pulses, and one can see the 
process of their generation even inside the soliton. We 
also note that the single cycle pulses move faster than 
the electron ridge. 

The frequency of the fields in the reflected single-cycle 
e. m. pulses is up-shifted and their longitudinal size is 
much smaller than the size of the soliton, as it can be 
clearly seen from Fig. [3J These simulation results confirm 
the physical mechanism investigated analytically in Ref. 
0, as summarized in Sec. Ill Al Since the soliton is 
not exactly positioned at the crossing of laser pulse axes, 
the reflected pulse is not exactly directed along the x- 
axis. This is a consequence of the parabolic profile of the 
wakefield: as the pulse is reflected by the upper wing of 
parabola it propagates at an angle with respect to the 
x-axis. 

Let us now estimate the parameters of the reflected 
single-cycle pulse according to the results of the ID ana- 
lytic calculations, using the initial conditions of the nu- 
merical simulations. In this case 7^ « 2.63. The in- 
tensity of the soliton is Is = {as^s l^df 'h ~ 0.0161%, 
where as ~ 0.5 is the soliton amplitude. The reflected 
pulse intensity according to the ID analytic prediction 
is I^ f = j^ h I s ss 0.28/i. The reflected pulse ampli- 
tude, as seen from the results of numerical simulations, is 
a re f ra 0.14. Then I ref = (2 1 2 h a ref n s /uJ d ) 2 h « 0.23Ji, 
which is in reasonable agreement with the ID theoretical 
analysis and reproduce well the predicted scaling. 



III. REFLECTION OF THE 
ELECTROMAGNETIC FIELD OF THE 
ELECTRON VORTEX 

A. One-dimensional theory 

Ultrashort electromagnetic pulses can also be gener- 
ated in the interaction of a breaking wake plasma wave 
with different types of coherent structures besides soli- 
tons, such as electron vortices, wake fields, etc. In this 
section we consider the interaction of a wakefield with an 
electron vortex. For the sake of illustration, in the fol- 
lowing we will employ a simplified representation of their 
spatial structure. 

An electron vortex is characterized by a quasi-static 
magnetic field generated by the current of electrons cir- 
culating in a plasma with steady ions. Contrary to hy- 
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FIG. 4: The electric E re f and magnetic B re f fields of the 
reflected pulse in the vortex-wakefield interaction as obtained 
from the one-dimensional theory vs. the spatial coordinate x 
at t = 0, for j p h ~ 10. The fields are measured in units of the 
initial vortex field and x is measured in units of c/u} pe - 



drodynamic vortices in ideal fluids, electron vortices in a 
plasma have a characteristic spatial scale that is given 
approximately by the plasma collisionless skin depth 
d e = c/uj pe . The generation of electron vortices, together 
with their associated magnetic field, by ultraintense laser 
pulses in plasmas is observed routinely in computer sim- 
ulations [111 (see also HH3 and [H). 

For the sake of simplicity we consider a ID configu- 
ration, which corresponds to magnetic "ribbons" as dis- 
cussed e. g. in [2]], [2^, take the magnetic field to be 
perpendicular to the direction of the wake wave and as- 
sume its profile to be gaussian, i. e. 

B = B exp(-x 2 /(f e ). (5) 

Following the procedure proposed in Ref. , we perform 
the Lorentz transformation to the reference frame where 
the wakefield is at rest and the magnetic field of the elec- 
tron vortex appears as an incident e. m. pulse. Then, 
we Fourier transform the incident pulse and use the re- 
flection and transition coefficients derived in Ref. Q . By 
inverting the Fourier transform and by performing the 
Lorentz transformation back to the laboratory frame, we 
obtain that the shape of the reflected pulse is given by 
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(see Fig. 0) 
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where <p = 2 1 / 2 / (3 ph JJ* , i> = ±V2j% 2 uj pe (t 

X 

erf(x) = -4= J exp(—t 2 )dt is the error function. As ex- 
o 

pected, the field of the reflected pulse scales for large j p h 
as 7^ 2 , and thus its intensity scales as "i ph - However the 



reflected pulse width has a more complicated scaling. As 
can be clearly seen from Fig. where the electric and 
magnetic fields of the reflected pulse are shown, we have 
two different scales ^h ph and ^/j p [ 2 , and the front of 
the reflected pulse is much steeper than its tail. This 
behaviour can also be deduced from Eq. I© since the 
width of the front of the pulse is determined by the ex- 
pression in square brackets where the argument of the 
error function scales as \/j^ h . On the contrary, the form 
of the tail is determined by the factor exp (—ip) where tp 



is proportional to 7, 



.3/2 
ph ' 



B. Two-dimensional PIC simulation 

The interaction of a wake wave with an electron vor- 
tex, associated with a quasi-stationary transverse mag- 
netic field, is simulated in a box with size 50 Ad x 80 A^, 
including absorbing edges of thickness 2Ad- The re- 
sults are shown in Figs. 15171 The electron vortex is 
prepared by introducing a magnetic field B z which in- 
creases gradually from t = to t = 10. After t = 10, 
a self-consistent quasi-stationary electron fluid vortex is 
formed with its corresponding magnetic field distribu- 
tion B z = 0.066 exp ( - (x - 20) 2 /4 - y 2 /4), Fig. 
The shape of the plasma slab is the same as in the pre- 
vious simulation (Sec. Ill Bp except that its transverse 
size (along the y-axis) is 70Xd and the maximum electron 
density is n e = 0.04n cr . The number of quasi-particles 
is 3.4 x 10 6 . A Gaussian, linearly polarized (E in the 
z-direction) , driver laser pulse with amplitude = 1.2 
and FWHM size 2\d x 30Ad, starts from t = 10 from the 
left boundary, propagates along the x-axis, and induces 
a plasma wakeficld which interacts with the electron vor- 
tex. The phase velocity of the wakefield is v p h ~ 0.965, 
so that the Lorentz factor is 7 P /, w 3.81. A large waist of 
the driver laser pulse is chosen to make the interpretation 
of the reflection easier. 

Fig. Elshows the ^-component of magnetic field and the 
electron density during the interaction. The wake wave 
itself has a transverse magnetic field B' z which arises due 
to the wake wave curvature. The wake wave bending 
implies that, in addition to the longitudinal oscillatory 
motion, electrons move towards the wake wave axis, thus 
inducing a weak transverse magnetic field which is max- 
imum on the periphery, in the regions of the electron 
density compression. Eventually, this effect results in 
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FIG. 5: The electric field component E z (color-scale) and 
contours of the magnetic field component B z for values B z = 
1,2,3,4,5,6 x 10 -2 (a) and the electron density n e (b) at 
t — 27 x 2n/u. The dashed circle denotes the electron vortex 
location. 



a so-called transverse wave break of the wakefield |24| . 
Although the magnetic field of the wake wave does not 
prevent us from distinguishing the e. m. pulses reflected 
from the vortex easily, nevertheless, in order to improve 
the presentation, in Fig. 0we show the quantity B z — B' z , 
where B z and B' z are the magnetic field components seen, 
respectively, in two simulations - one, as described in this 
section, and the other performed with the same parame- 
ters but without the vortex (i. e., with no initial magnetic 
field). Thus a "pure" reflection of the e. m. field of the 
vortex, resulting in formation of single-cycle pulses, can 
be seen. In Fig. 0the magnetic field component B z and 
electron charge density n e are shown along the x-axis 
at y — 0, where the magnetic field of the wake wave is 
exactly zero. 

As in the case of the interaction of the wake wave and 
a soliton, the electron density cusp is substantially dis- 
torted as it moves through the vortex and it recovers 
after passing the vortex, Figs. H3 El As in the soliton 
case this transient distortion of the cusp does not affect 
the shape of the resulting single-cycle pulses. We also 
note that the electron vortex is distorted due to modula- 
tions of the electron density in the wake wave. The size of 
the vortex and the magnitude of corresponding magnetic 
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FIG. 6: Interaction of the wake wave with the electron vortex, 
t — 47 x 2h/uj. (a) The difference between the magnetic 
field component B z and B' z , where B' z is obtained from the 
simulation without the vortex, (b) The electron density. 
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FIG. 7: Reflection of the electromagnetic field of the electron 
vortex by electron density cusps. The magnetic field compo- 
nent B z and electron density n e at t — 47,50 x 2n/u> along 
the :r-axis at y = 0. The first and second reflected pulses are 
marked with "1" and "2" , respectively. 



field change so that the angular momentum and vorticity 
of the electron fluid are preserved. 



IV. REFLECTION OF THE 
ELECTROMAGNETIC FIELD OF THE PLASMA 
WAKE WAVE 

A. One-dimensional theory 

Let us now consider the interaction of two wake plasma 
waves which arc oriented perpendicularly with respect to 
each other. We assume that the amplitude of one of 
the two waves (the first) is much smaller than that of 
the other (the second) and thus neglect the action of 
the first on the second wave. In this scheme the first 
wake wave provides the electric field which is (partially) 
reflected by the relativistic mirrors represented by the 
electron density spikes of the second wake wave. 

We assume that the electric field of the first wake wave 
is directed along the j/-axis and its spatial dependence 
along the x-axis, i. e. along the propagation direction 
of the electron density spikes of the second wave (the 
mirrors), can be modeled as a two-step function with a 
constant amplitude plateau in between. We can choose 
the parameters of the e. m. pulse that generates the 
first wake wave in such a way that the wavelength of the 
latter at breaking is much larger than its transverse size 
(the size of the amplitude plateau) which is assumed to 
be given by the characteristic scale of nonlinear plasma 
structures d e = c/uj pei i. e., such that 2^J2^ p hc/iTui pe 3> 
d e = c/uj pe . In this case the electric field of the first wave 
during its reflection from the second wake wave can be 
taken to be essentially a function of x only and to be of 
the form 



E w = E 9(x)9(c/uj pe - x). 



(7) 



Following the same procedure as in the previous section, 
we obtain for the electric field of the reflected pulse 

Eref = 2j ph E a 6(t - x/c) [1 - exp(-V0 



1 — exp 




(t-x/c) 



2uj 



(see Fig. [SJ. Here we also see the interplay of two scales: 

l/j ph and l/7^ 2 as in Fig. 21 The first one is connected 
with the width of the reflected pulse front and follows 
from the distance between the arguments of the two thcta 
functions. The second one determines the tail width and 
arises from the exp(— tf>) term in Eq. JSJ. The reflected 
pulse field strength reaches its maximum value 



1 — exp 




(9) 
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first scale is a consequence of the Doppler effect and is 
connected with the frequency and wavenumber upshift. 
It manifests itself in the frequency of the reflected pulse 
from a soliton and in the scale of the front parts of the 
reflected pulses from a vortex and a wakefield. The sec- 
ond scale manifests itself, in the case of the soliton, in 
the size of the reflected pulse envelope and, in the case 
of the vortex and of the wakefield, in the decay length of 
the reflected pulse tail. This second scale is due to the 
preferential reflection of the low frequency part of the 
field, i. e., to the fact that the reflectivity of the mirror 
increases and tends to unity for e. m. radiation with frc- 

1 /2 

quency of the order of j p ' h ui pe in the proper frame of the 

mirror (or 7p/( 2 w pe , in the laboratory frame). 

Finally, we should also note that ultrashort e. m. 
pulses can be generated when a wake wave interacts with 
plasma in a self-focusing channel |25|. see also, e. g., j2|| 
and references therein. In a self-focusing channel, an elec- 
tric field is present due to charge separation and can be 
transformed into an ultrashort e. m. pulse, similarly to 
the soliton, vortex or wake wave. 



-150 1 1 1 1 1 1 1 1 1 
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x 

FIG. 8: The electric E re f and magnetic B re f fields of the 
reflected pulse in the wakefield-wakefield interaction vs. the 
spatial coordinate x at t — for "/ p h = 10. The fields are mea- 
sured in units of the initial field strength of the first wakefield; 
x is measured in units of c/Lu pe . 



at t — x/c — 1 / (2jp h uj pe ) . For large values of j p h we 
obtain: 

E max = 4V2 7 ^ 2 £ - (10) 

Thus the amplitude of the e. m. field scales as j ph , 
similarly to the case of the reflection of the soliton and of 
the electron vortex from a breaking wake plasma wave. 
This scaling can be explained as follows. First, the e. m. 
field after reflection by the moving mirror acquires a fac- 
tor proportional to ^ h due to Doppler effect. This can 
be easily shown with the two Lorentz transformations - 
to the proper reference frame of the mirror and back, to 
the laboratory reference frame (the vector potential rep- 
resenting the e. m. field is perpendicular to the direction 
of propagation of the mirror and thus remains unchanged 
while the frequencies and wavenumbers are upshifted by 
a factor 7^). Second, the reflection coefficient of the 

— 1/2 

electron density spike is proportional to 7 A . These 

two facts result in the overall 7^ 2 factor. Actually, as 
discussed in the description of Figs. 0] and 00 the spatial 
scaling of the reflected pulse is not simply accounted for 
by the frequency and wavenumber upshift. In these two 
figures we observe two scales: 1/7^ and l/7p{ 2 - The 



B. Two-dimensional PIC simulation 

In order to show the interaction of a wake wave with 
another (weak) wake wave, we performed a simulation 
with the same simulation box and plasma distribution as 
in the previous cases (Sec. IIII Bp except for the maximum 
electron density which is now n e — 0.01n cr . The results 
are shown in Figs. 1911 II Two laser pulses, which are 
linearly polarized (E along the z-axis) and propagate in 
perpendicular directions, induce two plasma wakes, Fig. 
The first, auxiliary, laser pulse with Gaussian shape, 
amplitude a a = 0.7, wavelength A a = Xd and FWHM 
size (length x waist) 2.5Aci x lOA^, starts at t = from 
the bottom boundary, propagates in the direction of the 
y-axis at x = 20 and induces a weak wakefield with a 
longitudinal electric field directed along the y-axis. In 
Fig. [5] this weak wake wave is seen as a vertical periodic 
structure. The second laser pulse (driver) is switched on 
at t = 70 on the left boundary. Its amplitude is = 2 
and size is 2 Ad x 50Xd- The driver pulse has a Gaussian 
profile in the direction of the x axis; along the y-axis 
its shape is a smooth centrosymmetric, piecewise poly- 
nomial with sizes of adjacent parts 8,34,8. The central 
part is constant; the profile of the first slope is defined by 
(3 — 2r;)r; 2 , 77 = y/8. We chose such the shape of the driver 
laser pulse so as to ensure the transverse homogeneity 
of the induced wakefield in the region —17 < y < 17, 
and thus to obtain a plasma wake with zero transverse 
electric field in this region. In other words, the quasi- 
one-dimensional part of the driver induces the wakefield 
which is also quasi-one-dimensional. With such a config- 
uration, in the distribution of the E y component inside 
the region —17 < y < 17, we see, in principle, only the 
electric field of the weak wakefield, induced by the auxil- 
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FIG. 9: (a) The electric field component E y (color-scale) and 
contours of the component E z for values E z — —0.5 (dashed) 
and +0.5m e wc/e (solid) and (b) the electron density n e at t — 
85 x 2-k/uj. The weak wake wave, seen as a vertical periodic 
structure, is generated by the auxiliary laser pulse which has 
already left the simulation box. In the wake from the driver 
laser pulse, the E y is zero in the region —10 < y < 10, which 
corresponds to the homogeneous part of the driver. 

iary laser pulse. The wakefield from the driver has phase 
velocity v p h ~ 0.982 and the Lorentz factor j p h ~ 5.32. 

Fig. show the y-component of electric field and 
electron density at three moments of time with period 
approximately equal to the period of the plasma wake 
wave induced by the driver laser pulse. The driver is 
represented by contours of the ^-component of electric 
field. In Fig. ^] the electric field component E y and 
electron charge density n e are shown along the x-axis at 
y = 0. 

Again we note that, as in the case of the interaction 
of the wake wave with a soliton and vortex, the electron 
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FIG. 10: Interaction of the wake wave with the weak wave 
wave, shown in the portion of the simulation box, —10 < 
y < 10, corresponding to the quasi-one-dimensional part of 
the wakefield from the driver. The electric field component 
E y (a,c,e) and electron density n e (b,d,f) at t = 100 (a,b), 
110 (c,d) and 120 x 2ir/ui (e,f). The driver laser pulse is 
represented by contours of the electric field component E z in 
the same way as in the previous figure. 



density cusp is distorted inside the weak wake wave, but 
is almost restored when it moves outside, Figs. ^3 ITT1 
Even though the electron density cusp seems curved when 
it is under the influence of the weak wake wave, the single- 
cycle pulse reflected by this cusp appears to be flat. 
As in the cases of other nonlinear coherent structures 
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FIG. 11: Reflection of the electromagnetic field of the weak 
wake wave by electron density cusps of the wake wave from the 
driver. The electric field component E y and electron density 
n e at t = 115, 120 x 2n/u) along the :r-axis at y = 0. The first 
and second reflected pulses are marked with "1" and "2", 
respectively. 



described above, each density cusp reflects single-cycle 
c. m. pulse whose frequency is up-shifted and whose in- 
tensity is increased due to Doppler effect. A complica- 
tion, which is not included in the presented ID theory, 
arises from the fact that the weak wake wave, which is 
(partially) reflected by the "mirrors" of the wake from 
the driver, has almost the same phase velocity as these 
mirrors. 



presented here confirm and extend the analytical results 
obtained in Ref. 0]. 

The modulations of electron density in the strong wake 
wave, which is close to the wave-breaking regime, have 
the shape of spikes. Each spike acts as a semi-transparent 
mirror moving with a relativistic velocity, corresponding 
to the phase velocity of the wake wave. Such a mirror 
partially reflects the electromagnetic field of a coherent 
nonlinear structure and thus generates an electromag- 
netic pulse. As predicted in Ref. Q, the reflected pulse 
consists of a single cycle oscillation and, as compared to 
the e. m. field of the coherent nonlinear structure, the 
reflected pulse has an up-shifted frequency and an in- 
creased intensity. The reflected pulse intensity occurs as 
a result of frequency up-shift, due to Doppler effect, and 
because of the parabolic profile of the wake wave. 

Using an analytical approach, we have shown that in 
the three cases considered here of coherent nonlinear 
structures - a sub-cycle soliton, an electron vortex and a 
wake wave - the amplitude of the e. m. pulse, reflected 
by a relativistic flying mirror, scales as 7p{ , due to a 
similarity of the reflection process in all three cases. The 
reflection leads to a frequency and wavenumber upshift 
which scales as jz h and to the formation of an additional 

spatial scale proportional to l/7p^ 2 - 

Since in all the above cases of wakefield interaction 
with nonlinear coherent structures in a plasma, single- 
cycle e. m. pulses are emitted with a characteristic fre- 
quency, duration and polarization, their emission repre- 
sents an important process to be used for diagnostics of 
laser plasma interactions. 



CONCLUSION 



Our analytical model and two-dimensional particle-in- 
cell simulations show that during the interaction of coher- 
ent nonlinear structures (such as sub-cycle solitons, elec- 
tron vortices and wake Langmuir waves) with a strong 
wake wave in a collisionless plasma a train of single-cycle 
intense electromagnetic pulses is generated. This effect 
can be exploited in order to produce ultra-short intense 
e. m. pulses with presently available lasers. The results 
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